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MiSTRACT To clarify phylogenetic relationships among the decabraehia cephalopods, especially the family Sepiidae and Sepiol- 
idae, mitochondrial cytochrome c oxidase subunit 1 (COl) gene and partial 16S rRNA gene were sequenced tor 13 species. Phylogenetic 
analyses were performed by the distance and parsimony method, Coleoids were divided into 2 main lineages, Octobrachia and 
Decabraehia (including Sepiida. Sepiolida, and Teuthida). In all phylogenetic trees, the monophyly of the Sepiolidae and Sepiidae was 
supported well, but their rank and position within the Decabraehia were not clear. Based on partial COI rDNA and its amino acid 
sequences. Parsimony analyses showed Sepiolidae. Sepiidae, and the ( Loligo chinensis , Chtenoteryx sicnla ) cladc from Teuthida were 
in the same level. Compared with Sepiidae, Sepiolidae was more closely related to Teuthids using I6S rDNA sequences. We are 
inclined to support the current classification: Sepiolidae and Sepiidae belong to different Orders. According to the phylogenetic 
analysis, the 2 genera ( Sepiella and Sepia) from the Sepiidae can be distinguished (7X% neighbor-joining (NJ); 64% maximum 
parsimony (MP) in I6S rRNA gene), but do not have visible boundline using COI gene and its amino acid data. This suggests that COI 
gene may be much fitter to analyze cuttlefish phylogeny at a high taxonomic level (i.e.. family), and 16S rRNA gene could be used 
as a precious tool to analyze taxonomic relationships ai the genus level. 
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INTRODUCTION 

According to Voss (1977), the Sepioidea includes Spirulidae, 
Sepiidae, Sepiolidae, ldiosepiidae, and Sepiadariidae. To date, the 
taxon system is still used in China. The position of the Sepiolidae 
as a sister group of the Sepiidae and the monophyly of the squids 
are, however, questioned by some scholars (e.g.. Berthold & En- 
geser 1987, Clarke 1988, Boletzky 1999). 

Molecular information from DNA data has been used to clarify 
the relationships among Cephalopods since the 1990s. The nucle¬ 
otide divergence between sequences of one gene or a portion of 
gene from the mitochondrial or nuclear genome can be analyzed 
phylogenetically. Previous DNA sequence diversity and phyloge¬ 
netic relationships of octopods have been investigated using the 
mitochondrial cytochrome c oxidase subunit 1 (mtCO I), mtCO II, 
mtCO 111. 16S rRNA gene (Carlini & Graves!999, Carlini et ai. 
2001, Bonnaud ct ai. 1996, Bonnaud et al.1997, Soller et al. 2000, 
Allcock & Piertney 2002, Piertney et al. 2003). Anderson (2000) 
sequenced 2 mitochondrial genes {16S rRNA and CO I) to clarify 
loliginid phylogeny. Sequence analyses from the 3' end of the mt 
1-rRNA (I6S) gene of decapod cephalopods have shown that this 
portion of gene was a useful tool for taxonomic relationships at the 
infrafamilial level (Bonnaud et al. 1994). The COI gene for phy¬ 
logenetic analysis of the coleoid cephalopods exhibited a high 
degree of nucleotide sequence variability, with one half of the sites 
varying in at least one taxon; COI amino acid sequences were 
highly conserved, but were useful in determining basal-level rela¬ 
tionships among the Coleoidea (Carlini et al. 2001). 

The cephalopods, especially decapods, are an important and 
valuable fishery resource in China. South Korea, and Japan (Nesis 
& Kir 1982, Okutani 1995). For example. Sepia esculenia , Sepiella 
maimlroni and Loligo chinensis are all high commercial species 
(Dong 1991). A great deal of fundamental and applied researches 
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on cephalopods, such as fauna, systematics. morphology, embry¬ 
ology, population genetics, and biodiversity have been completed 
since 1960s (e.g.. Lee 1963, Lee 1983, Dong 1993. Lu 1998. Lu 
2000, Zheng et al. 2001a, Zheng et al. 2001b, Zheng et al. 2004). 
There are, however, few documents referring to the molecular 
evolution and phylogenetics of cephalopods living in the coastal 
waters of China. In this study 13 cephalopod species are analyzed 
by the sequence comparison of 16S rRNA and COI gene. The 
phylogenetic trees are reconstructed. The taxonomic relationships 
of decabraehia species are discussed. 

MATERIALS AND METHODS 

Taxon Selection 

Eight decabraehia species and two octopus species were 
sampled for 16S rRNA gene sequence analysis. Representatives of 
three other cephalopods {Euprynuui scolopes, Bonnaud et al. 1994: 
Chtenopieryx sicnla , Grimpoteuthis sp., Anderson 2000) were 
added to cladistic analysis. 

Eight decabraehia species and one octopus species were 
sampled for the COI gene sequence analysis. Representatives of 4 
other cephalopods ( Nautilius pompilins , Carlini & Gravesl999; 
Enprynma scolopes , Bonnaud et al. 1994, Chtenoptery.x sicnla . 
Grimpoteuthis sp. Anderson 2000) were added to cladistic analy¬ 
sis. 

Details on the taxonomic position (following Voss 1977) and 
origin of the 13 species studied are presented in Table I. 

DNA Extraction , Amplification , and Sequencing 

Total genomic DNA was extracted using a CTAB method 
modified from Winnepenninckx et al. (1993). Regions of the 16S 
and COI genes were amplified by PCR. The 16S primers se¬ 
quences were D16SAR 5'-CGC CTG TTT AHY AAA AAC AT 
3', DI6SBR 5'-CCG GTC TGA ACT CAG MTC AYG 
T-3'(Anderson 2000). The primers used for the amplification of 
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TABLE I. 

Classification (Voss 1977) of cephalopod taxa included in this study. 


Classification 

a References 

Locality 

GcneCOI/16S# 

Phylum MOLLUSCA 

Class CEPH ALOPDA 

Subclass NAUTILOIDEA 

Family Nautilidae 

Nautilius pompilius Linne. 1758 

CA 

Waikiki Aquarium 

AF000054/— 

Subclass COLEOIDEA 

Order SEPIOIDEA 

Family Sepiidae 

Sepia aculeata Orbigny, 1848 

ZH 

East China Sea 

AF350494/AF369113 

S. esculenta Hoyle, 1885 

ZH 

Yellow Sea of China 

AF359554/AF369114 

5. latimanus Quoy & Gaimard, 1832 

ZH 

South China Sea 

AYI85506/AF369116 

S. pharaonis Ehrenberg, 1831 

ZH 

South China Sea 

AF359555/AF369117 

S. rohsoni (Massy, 1927) 

ZH 

South China Sea 

AF350495/AF369957 

Sepiella maindroni de Rochebrune. 18S4 

ZH 

East China Sea 

AF340032/AF369118 

Family Sepiolidae 

Euprymna berryi Sasaki, 1929 

ZH 

South China Sea 

AF350493/AF36911 ft 

E. scolopes Berry, 1913 

AN/BO 

Hawaii, USA 

AF075417/X79592 

Order TEUTHOJDEA 

Suborder MYOPSIDA 

Family Loliginidae 

Loliga chinensis Gray, 1849 

ZH 

South China Sea 

AY 1S5505/AF369955 

Suborder OEGOPSJDA 

Family Chtenopterygidae 

Chtenopteryx siciila (Verany, 1851) 

AN 

Pacific Ocean 

AF075416/AF110097 

Order OCTOPODA 

Suborder ClRRATA 

Family Cirroteuthidae 

Grimpotcuthis sp. 

AN 

Monterey Bay, CA 

USA 

AF075419/AF110100 

Suborder 1NCIRRATA 

Family Octopodidae 

Octopus ocellatus Gra>. 1849 

ZH 

East China Sea 

AF346S54/AF369111 

O . variabilis (Sasaki. 1929) 

ZH 

East China Sea 

—/AF369112 


a References: AN, Anderson F. E. (2000); CA, Carlini, D. B. and Graves. J. E. (1999); BO, Bonnaud, L. Boucher-Rodoni R, Monnerot M. (1994). ZH, 
the authors. #Genbank accession numbers. 


partial COl gene were: HC02198 (5'-TAA ACT TGA GGG TGA 
CCA AAA AAT-3') and LCOI49I (5'-GGT CAA CAA ATC 
ATA AAG ATA TTG-3') from Folmer et al (1994). Amplifica¬ 
tions were performed under the following conditions: 120 s at 
94 °C, then 30 or 35 cycles (COl and 16S, respectively) of 40 s at 
94 °C, 1 min at 50 °C, and 1 min at 72 C. A total volume of 25 
p.L reactions consists of 0.5 units of Taq (TaKaRa), 0.5 pM each 
primer, 0.2 pM each dNTP. 2.5-3.5 mM MgCU, 2.5 pi of lOx 
buffer supplied with Taq and 4 pi (30-50 ng) of template DNA. 

The amplified fragments were purified by the PCR fragment 
recovery Kit (TaKaRa). Purified products were sequenced directly 
using the AB1 PRISM BigDye Terminator cycle sequencing Ready 
Reaction Kit and AmplicTaq DNA polymerase with AB1 PRISM 
377X1. DNA sequencer (Applied Biosystem Inc.). 

Phylogenetic .4 nalyses 

All of the initial sequences were aligned by ClustalX v 1.8 
(Thompson et al. 1997). Amino acid sequences of the COl gene 
were translated using the GENEDOC (Nicholas et al. 1997). The 
mean nucleotide composition, proportion of transition (ts) by num¬ 
ber of the total base substitutions (ts + tv, transitions + transver¬ 


sions) were calculated in MEGA 2 (Kumar et al. 2001). The phy¬ 
logenetic analyses were performed on both aligned nucleotide se¬ 
quences and on amino acid sequences, using distance and 
parsimony methods included in the packages MEGA 2 (Kumar et 
al. 2001) and PAUP4.0b6 (Swofford 2000). Statistical confidence 
of a particular cluster of sequences w'as evaluated by the bootstrap 
procedure (1.000 resampling replicates). 

RESULTS 

Genetic Variation of Partial COl Gene Sequence 

The percentage proportion of 4 kinds of base pairs (A, C, G. T) 
was compared with each other. There was no remarkable differ¬ 
ence betw een these species examined. The average content of base 
pairs in the nine species examined was 28.12%. 17.51% , 15.71%. 
and 38.60%. respectively. A + T content was also up to 66%, 
which is in the same range as in the reference species. High A + 
T content was an obvious characteristic of cephalopod mtDNA 
sequences. No gaps were found in all of the sequences we ana¬ 
lyzed. A total of 281 nucleotide positions were found to be variable 
(42.4%). One-hundred and ninety-four variable sites were found at 
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the third position of codon triplet and the percentage was approxi¬ 
mately 69% of all variable nucleotides. 

The genetic distance ranged from 0.000 to 0.202 among the 
Sepiidae. Sepiella maindroni could not be separated from other 
species of the family using distance method, although it belonged 
to Sepiella in the traditional taxonomy. The closest taxon with 
respect to genetic distance was Sepia latimanus (0.000), and the 
most distant was Sepia aculeata (0.202). Compared with the data 
between Sepiidae and Sepiolidae, the distance was 0.191 ± 0.008 
(mean ± SD), slightly closer than between Sepiidae and Teuthida 
(mainly mentioned Loliginidae and Chtenopterygidae) (0.206 ± 
0.007) as well as between Sepiolidae and Teuthida (0.215 ± 
0.013). Nautilus pompilius , belonging to Nautilida, was the most 
distant species from the other 12 cephalopods analyzed (0.288 
-0.349). 

According to the amino acid sequences, the genetic distance 
was in the range of 0.000 to 0.0056 within the Sepiidae. Sepiella 
maindroni was grouped with the other cuttlefishes, even though 
they belong to a different genus in the traditional taxonomy. The 
larthe.st distance between Nautilus pompilius and other cephalo¬ 
pods (0.119 -0.178) indicated that the relationship between them 
was the most remote. 

Genetic Variation of Partial I6S rRNA Gene 

The length of 16S rRNA gene sequences of the 10 species 
examined ranged from 500 to 514 base pairs (average 507.5 bp). 
The average of percentage proportion of A, C, G, and T was 
34.25%, 9.83%, 17.93%, and 37.99%, respectively. The content of 
A + T is up to 72%. The result is similar to the above research of 
the COl gene sequences for cephalopods. 

In the Sepiidae, the genetic distance of 5 species of Sepia 
ranged Irom 0.030 to 0.045, and become a bit higher with Sepiella 
(0.059 -0.077). The distance between Sepiidae and Sepiolidae 
(0.085 ± 0.009) was not significant different from that between 
Sepiidae and Teuthida (0.090 ± 0.013) or Sepiolidae and Teuthida 
(0.089 ± 0.012). The longest distance (0.173) was observed be¬ 
tween Sepiolidae (Euprynma scolopes) and Octopoda (Octopus 
variabilis). 


Phytogeny of Cephalopods Rased on COl and I6S rRNA 
Gene Sequence 

Phylogenetic Analysis Based on Partial COl and Amino 
Acid Sequences 

The topology structure of phylogenetic trees was similar using 
the COl and their amino acid sequences, whatever NJ trees and MP 
trees. The monophyly of the Sepiolidae (98% NJ, 78%MP in Fig. 
1; 90% NJ, 74% MP in Fig. 2) and Sepiidae (50% NJ, 58% MP in 
Fig. 1: 93% NJ, 87% MP in Fig. 2) was well supported. According 
to NJ analyses. Coleoids were divided into 2 main lineages, Oc- 
tobrachia and Decabrachia (including Sepiidae, Sepiolidae, and 
Teuthida). In the Decabrachia clade, 6 species from the Family 
Sepiidae grouped with 2 species from the Order Teuthida at first, 
then grouped with the Family Sepiolidae (NJ trees in Fig. 1 and 2), 
though low levels of bootstrap support were obtained (see asterisk 
place in Fig. 1 and 2). Parsimony analyses indicated that the De¬ 
cabrachia was divided into 3 parallel clades (88% MP in Fig. 1: 
96% MP in Fig. 2). Within the Sepiidae clade, Sepiella ( S . main¬ 
droni) and Sepia could not be divided clearly; and the groups from 
5 species in Sepia were not consistent with morphologic evidence, 
which belong to 3 different species complexes (Khromov 1998). 

Phylogenetic Analysis Based on Partial 16S rRNA Gene Sequence 

Both MP and NJ analyses provided strong bootstrap support 
values for the monophyly of the Coleoid cephalopods (98% MP 
and 98% NJ), and the monophyly of the Family Sepiidae (64% MP 
and 78% NJ) as well as Family Sepiolidae (95% MP and 96% NJ) 
(Fig. 3). Within the Decabrachia clade, Sepiolidae was more 
closely related with Teuthida (54% MP and 72% NJ). Sepiella and 
Sepia were clearly separated (64% MP and 76% NJ). Though the 
topology structure (NJ tree) in Sepia clade was consistent with the 
classilication of species complexes of Sepiidae, bootstrap support 
values were very low (<50%). 

DISCUSSION 

Previous investigations of coleoid systematics have attempted 
to determine relationships within the Sepiida, Octopoda. and 
Teuthida through phylogenetic analysis of morphologic and mo¬ 
lecular character data. Most recently, Carlini and Graves (1999) 
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Figure t. NJ phylogenetic tree and MP tree based on COl gene data. Boldfaced numbers above branches are bootstrap support values (1,000 
replicates). Asterisk indicated bootstrap values less than 50%. Nautilus pompilius was used as distant outgroup species. (1) Sepiidae. (2) Teuthida, 
(3) Sepiolidae, (4| Octopodida. 
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Figure 2. NJ phylogenetic tree and MP tree based on amino acid of COI gene data. Boldfaced numbers ahove branches are bootstrap support 
values (1,(100 replicates). Asterisk indicated bootstrap values less than 509c. Xautiluspompilius was used as distant outgroup species. (1 ) Sepiidae, 
(2) Teuthida, (3) Sepiolidae, (4) Octopodida. 


conducted a higher-level analysis of the coleoid cephalopods. Con¬ 
sistent with them, our results suggest that the coleoids can be 
divided into 2 main lineages: the Octobrachia and the Decabrachia. 

According to the classification of cephalopod taxa (Voss 1977). 
the Order Sepioidea was comprised of 5 families Spiru]idae, Se¬ 
piolidae. Sepiidae, Sepiadariidae. and Idiosepiidae (Fig. 4). Al¬ 
though the COI and 16S rDNA data well supported the monophyly 
of the Decabrachia, we could not confirm the validity of the order 
Sepioidea. Several studies have questioned this order (Clarke 
1988, Bonnaud et al. 1996. Bonnaud et al 1997. Sweeney & 
Roper 1998, Car)ini & Graves 1999) (see Fig. 4). Bonnaud et al 
(1994) analyzed phylogenv of decapod cephalopods based on 3' 
end of 16S rDNA nucleotide sequences, and demonstrated that 
Sepiolids were clearly excluded from the order, and the position of 
the Spirulidae required further clarification. Referring to paleon¬ 
tology and neontology evidences of cephalopods, Clarke (1988) 
pointed out that the fiv e families belonged to two orders, that is the 
Order Sepiida was composed of Sepiidae. Spirulidae and Sepia¬ 


dariidae; Sepiolidae and Idiosepiidae belonged to the Order Sepi- 
olida. Sweeney and Roper (1998) listed the currently accepted 
classification of the recent Cephalopoda and addressed that the 
five families belonged to three orders. Boletzky (1999) indicated 
that they should be subject to four orders, and the Superorder 
Decabrachia included 5 orders. Spirulida, Sepiida, Sepiolida, Idi- 
osepiida, and Teuthida (see Fig. 4). Carlini and Graves (1999) 
pointed out that the Sepioidea were polyphyletic. In all phyloge¬ 
netic trees, the monophyly of the Sepiolidae and Sepiidae was 
supported well, but their rank and position within the Decabrachia 
were not clear. Based on partial COI and its amino acid sequences. 
Parsimony analyses showed Sepiolidae. Sepiidae as well as the 
{Loligo chinensis . Chtenoteryx siatla) clade from Teuthida were in 
the same level (S89F in Fig. I and 96# in Fig. 2). Compared with 
Sepiidae, Sepiolidae was more closely related to Teuthids (72# NJ 
and 54# MP) using I6S rDNA. We are inclined to support the 
current classification: Sepiolidae and Sepiidae belong to different 
orders. As for the rank and position of other families (e.g., ldi- 
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Figure 3. NJ phylogenetic tree and MP tree based on based on 16S rRNA gene data. Boldfaced numbers above branches are bootstrap support 
values (1000 replicates). Grimpoteuthis sp. was used as distant outgroup species. (1) Sepiidae, (2) Sepiolidae, (3) Teuthida, (4) Octopodidae. 
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Figure 4. Sepiid relationships with other ta\n according to various 
authors. Ordinal rank indicated by capital letters. 


osepiidae, Spirulidae. etc), such researches will he carried out 
further and testify these families' phylogenetic relationship. 

According to the COl gene and amino acid sequence data, the 
family-level relationship was exhibited clearly, but the genus 
Sepiella (5. maindroni )was grouped with other species from Sepia 
(see Fig. 1 and 2). Sepiella was obviously separated from the Sepia 
groupings based on the 16S rDNA sequences (see Fig. 3), which 
demonstrated that they belonged to two separate genera. It shows 
possibility that mtDNA sequence fragments (e.g., 16S rDNA. 
COl) with different evolutionary rate are fit for different taxa 
phylogeny. The 16S rDNA of Cephalopods is a precious tool to 
analyze taxonomic relationships at the genus level, and COl gene 
is fitter at the family level. Of course, 3 species complexes within 
the genus Sepia were not well supported in all the MP trees. The 
recognition of subgeneric rank in Sepia is not used because more 
analyses on morphologic, behavioral, biochemical, and molecular 
data are needed to resolve the question (Lu 2004, pers. comm.). 
We will examine more extensive sequence data of the species of 
Sepia for more refined information of species complexes or sub¬ 
generic rank. 
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